Topological transitions of gapless paired states in mixed-geometry lattices 
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Fermion pairing may coexist with magnetism in 
unconventional superconductors and nuclear matter, 
stabilized by non-BCS mechanisms. Remarkable sce- 
narios of exotic pairing with Fermi surfaces mismatched 
by chemical potential or mass difference have been 
suggested [HE], inspiring experiments in solid-state 
materials [5j [6] and ultracold Fermi gases [7HT0] . Here 
we propose fermionic species to be confined in lattices 
of different geometry and show how such distortion 
of symmetry leads to unconventional quantum states. 
In particular, we consider pairing of one species in 
a graphene-like lattice with another in a triangular 
sublattice. We find a rich phase diagram of multiband 
pairing with gapped and gapless excitations at zero 
temperature, implying that gapless superfluids can be 
stabilized in these systems. Quantum phase transitions 
between topologically distinct phases are predicted. 
Realization of mixed-geometry systems, e.g. using 
spin-dependent optical lattices [TTHT5] , provides a new 
dimension to the ultracold quantum simulator for 
exploring unconventional quantum states of matter. 

The standard BCS theory considers a perfectly sym- 
metric pairing in up- and down-spin particles with op- 
posite momenta lying on identical noninteracting Fermi 
surfaces. One of the fundamental issues in condensed 
matter physics is what would happen if this setting was 
broken, for instance, with mismatched Fermi surfaces 
between the two spin components. The Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state suggests Cooper pairs 
carrying finite center-of-mass momentum, causing a spa- 
tially oscillating order parameter [TJ[2]. The Sarma or 
breached-pair (BP) state describes polarized superfluid 
of BCS-type zero-momentum pairs separated from un- 
paired excess particles forming a Fermi sea [3l [4] . These 
scenarios remain elusive but inspiring as possible uncon- 
ventional pairing mechanisms. 

Ultracold Fermi gases are potentially an ideal testbed 
of quantum many-body physics because of their un- 
precedented controllability [16]. Loaded in optical lat- 
tices, fermionic atoms with tunable interparticle inter- 
action can construct the repulsive or attractive Hub- 
bard model [T7J [18]. One remarkable feature of ultra- 
cold gases is the tunability of asymmetry between the 
(pseudo-)spins associated with atomic internal states. 
Typical forms of spin asymmetry are population [7HTU] 
or mass [19 -21 imbalance between the spin components, 
or spin-dependent optical lattices realized for bosons [TT1 - 
[15] • The effects of anisotropic hopping or onsite poten- 
tials in a lattice [22j [23], as well as mixed-dimensional 



pairing in continuum [24, 25 , have been theoretically in- 
vestigated. Here we propose a fundamentally different 
type of asymmetry: a spin-dependent lattice where each 
spin component is separately loaded on a different un- 
derlying geometrical setup. 

We propose a mixed-geometry setup of a two com- 
ponent Fermi gas where up- and down-spin components 
are loaded on honeycomb and triangular lattices, respec- 
tively, as shown in Fig. []Jl. We argue that this is one 
of the simplest yet nontrivial mixed-geometry configura- 
tions. It is a superlattice of two spin-dependent sublat- 
tices A and B. Tuned on short-ranged interactions, this 
configuration can realize the Hubbard model with the on- 
site interaction selectively applied at A sites. The setup 
could be experimentally approached by generalizing the 
techniques in [TTHT5] or using two different atomic species 

jm. 

We consider the attractive Hubbard model for a 
fermion pairing problem in the mixed- geometry lattices. 
As described in Methods, the mean-field Hamiltonian can 
be written in momentum space as 
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which reveals a multiband system where up-spin parti- 
cles in band a = 1, 2 would pair with down-spin particles 
in band a = 3. The lattice effects of the mixed geom- 
etry are encoded in the noninteracting band dispersion 
£ a (k) sketched in Fig. [IJd and in the coupling gi^(k). 
With fixed interaction and hopping strengths, two con- 
trol parameters can govern the pairing in this system: 1) 
tuning chemical potentials can lead to the Fermi surface 
mismatch; 2) the onsite energy modulation e controls a 
relative strength of the couplings g\ and g<i. 

Figure [T]i shows the zero-temperature phase diagram 
for the case of the graphene-like up-spin bands (e = 0). 
We find that the diagram is divided into three main ar- 
eas indicating the normal, gapped, and gapless phases. 
While the normal phase is simply indicated by a vanish- 
ing order parameter A = 0, a paired state with nonzero 
A is discriminated by its single-particle excitation spec- 
trum to be in the gapped or gapless phases. The gapped 
phase is a fully paired state, similar to the BCS state. In 
the gapless phases, we find an exotic paired state which 
we call the incomplete-breached-pair (iBP) state because 
of a partial breach found in momentum distribution. 

Intriguingly, the iBP state can be further divided into 
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different subclasses by the topological arrangement of 
one or two Fermi surfaces, leading to Lifshitz transi- 
tions [26]. The iBP states undergo two types of tran- 
sitions where the Fermi surface is either vanishing (2-FS 
to 1-FS) or transforming to a topologically different sur- 
face (r-centered to K-centered). In this mixed-geometry 
system, while the former is a continuous transition, it 
turns out that the latter can be continuous or discontin- 
uous transitions which meet at a multicritical point (see 
Figs. [T^ij2^i). This end point of the discontinuous transi- 
tion may be explained as the marginal quantum critical 
point of a Lifshitz transition in a two-dimensional inter- 
acting system [27] . 

The nature of the iBP states can be characterized 
by momentum distributions n^^(k) and a momentum- 
resolved order parameter shown in Fig. [2] While 
the jumps in n^^(k) and A& indicate the presence of 
the Fermi surfaces, the breach surrounded by these dis- 
continuities shows clear contrast to that of the conven- 
tional BP state in momentum space. It turns out that 
the momentum-space breach in our iBP state shows a 
mixture of paired particles with finite A& coexisting with 
unpaired excess majority species, rather than a perfect 
phase separation suggested by the conventional BP state. 
However, the Luttinger's theorem still holds for a con- 
served quantity (here — n\) in the broken symmetry 
phase [28 , explaining the plateau of n^ — n^ = 1 observed 
between the Fermi surfaces. 

Controlled by chemical potentials or densities, topo- 
logically different iBP states arise with distinct Fermi 
surface topologies. It is not only that the number of the 
Fermi surfaces essentially changes between one and two 
as observed for instance at ^ = 0, but also the shape 
of the Fermi surface undergoes a topological transition. 
Figure [2] shows the topologically different Fermi surfaces 
found in this system and the phase diagram of the corre- 
sponding distinct iBP states. Note that the transitions in 
our phase diagram are driven by densities or chemical po- 
tentials, which is in contrast to the BEC-BCS crossover 
for imbalanced two-component gases where the transition 
into the BP state with one Fermi surface may occur in the 
BEC regime driven by interaction [29j [30] . In addition, 
compared with a noninteracting case where a topologi- 
cal change of the Fermi surface occurs simply at the van 
Hove singularity at |/i| = 1 (see Fig.JIj:), it is notable that 
the transition lines in the interacting system are very dif- 
ferent from \n\ = 1 as indicated by our diagram of the 
phases with T- and K-centered Fermi surfaces. 

The characteristics of the gapless phase can be also 
tuned by onsite energy modulation e. The effects of finite 
e are two- fold. First, it opens a band gap in the nonin- 
teracting up-spin bands. Second, it creates an imbalance 
between the interband couplings gi and #2 (see Eq. ([I])) 
as described in Methods. Consequently, for a large neg- 
ative e giving #2 ^> #1, the contribution of the up-spin 
upper band (a = 1) to the pairing can be significantly 



suppressed. Figure [3] displays the phase diagrams for 
e = —1.0,-1.4,-1.8, indicating that the gapless phase 
becomes less robust as |e| increases. This apparent con- 
nection between the coupling imbalance and the stability 
of the gapless phase implies that a well-balanced inter- 
band contribution may be an important factor for our 
mixed-geometry lattice in stabilizing the iBP phase. 

One of the main findings is the chemical-potential- 
driven emergence of the paired states with zero, one, 
and two Fermi surfaces. The importance of the multi- 
band contribution to the formation of these exotic paired 
states can be intuitively understood in the band struc- 
ture of the mixed-geometry system. Revisiting our mean- 
field procedures (see Methods), as illustrated in Fig. [4j 
there are three cases of zero, one, and two crossings al- 
lowed between the noninteracting particle-like and 
hole-like (—£4,) dispersions with varying chemical poten- 
tials. The interaction (#1,2) causes band repulsion and 
gaps at crossings, leading to the two different configu- 
rations of quasiparticle energy dispersions. Depending 
on the chemical potentials, the zero excitation level in- 
dicates the formation of zero, one, or two Fermi surfaces 
(see Fig. [4)3). 

The simple band picture can also clarify the arrange- 
ments of the distinct iBP states in the phase diagrams. 
For instance, considering the iBP state with two Fermi 
surfaces, Fig. indicates that there are two multiband 
configurations forming such states, the one in the lower 
band side (d) and the other in the upper band side (f). In 
contrast, the state with one Fermi surface occurs only in 
the upper band side (e). This classification in the band 
picture is consistent with the phase boundaries in the 
phase diagrams. In addition, the band picture indicates 
that the iBP states with two Fermi surfaces have a dif- 
ferent nature depending on whether they originate from 
the upper or lower band side. Indeed, for instance, the 
momentum distributions found at e = —1.8 show a fully 
polarized normal area of excess majority particles around 
K, which is in contrast to the fully paired area around K 
observed in its counterpart found at e = (see Fig. [3^) . 
Systematic study of finite temperature effects and quan- 
tum fluctuations is beyond the scope of this work. While 
the Fermi surface is strictly defined at zero temperature, 
we have tested that the main features are present at fi- 
nite temperatures within the mean-field theory, critical 
temperatures being typically of the order of a tenth of 
the hopping amplitude. 

Mixed-geometry lattice systems as we propose here ex- 
tend the realm of pairing and strong correlation physics. 
The attractive Hubbard model in the mixed honeycomb- 
triangular lattice that we have considered reveals a new 
exotic multiband paired state, the incomplete-breached- 
pair state, and shows a rich phase diagram of topolog- 
ically distinct phases. We have found that multiband 
pairing plays a key role in stabilizing the coexistence of 
pairing and magnetization, which may have an interest- 
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ing connection to multiband effects and magnetic fluctu- 
ations interplaying in solid-state superconductors such as 
iron-based materials and MgE>2. We argue that realizing 
mixed-geometry lattices would open a new dimension to 
the search for novel quantum states with ultracold gases 
as quantum simulators. 



METHODS 

Attractive Hubbard model. The Hubbard Hamil- 
tonian in the mixed-geometry lattice system considered 
becomes 

« = -*t E («it^t+ /i - c -)+etE^t+4E^t 
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where a) (a) and W (6) are fermionic creation (anni- 
hilation) operators in the sublattices A and B, respec- 
tively, and the density operators h a and h b denote a) a 
and tfb. The spin-dependent hopping occurs between 
nearest-neighboring sites in the up- spin honeycomb 
lattice jC^ and the down-spin triangular lattice while 
the hopping strengths and t± are set to be unity for 
simplicity. Introducing the control parameter e for on- 
site energy modulations, the onsite terms are chosen as 
= e/2 = — e/2, and = — 3 without loss of general- 
ity since other possible tunings of onsite energies can be 
absorbed into the spin-dependent chemical potentials ^ 
and /^i . The onsite interaction is chosen as U = 5. 
Mean-field theory. The mean-field order parameter 
is defined as A = U (a^a^). Neglecting fluctuations, the 
interaction term can be approximated as — U ^t^ti ~ 

-E,( A 4t4| + A*a^a it ) + \A\ 2 /U. Possible Hartree 
shifts can be incorporated into the chemical potentials 
and the onsite energy modulation e. Performing the 



Fourier transformation, = 



f ia where / 



is a or 6, the noninteracting part of the Hamiltonian 

(ai.bt.at) as 
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can be written in the basis of <£>t 
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— t^[cos fc^ + cos ^(fc^ + V^fcaj) + cos \{k y — V3k x ) — 3] — 
Hi. Accordingly, the order parameter can be rewrit- 
ten with the momentum resolution as A = A^ = 
U J^ki^-ki^k^) - We see that, in momentum space, the 



up-spin particles at different sites (A or B) are cou- 
pled. A more natural form of the noninteracting Hamil- 
tonian is obtained by choosing a new basis set that di- 
agonalizes The transformed fermionic operators 

k 

c^g and c a j: with the index a G {1,2,3} then represent 

the noninteracting energy bands where 1 and 2 

denote the upper and lower bands of the up-spin compo- 
nent, and 3 indicates the down-spin band. In the basis 
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the full mean-field Hamiltonian is 

A\ 2 /U 



(2) 



rewritten as n MF = EjJ + EjJ & (-*) 

with the matrix representation 

gt(k) 

S2(fe) -&(-*), 
where the diagonal terms are the noninteracting band 
dispersions £1,2 (&) = ±^/e 2 + |/^(/c)| 2 — /i^ and ^(fc) = 
and the interband couplings are derived as 



1 ± 



|/> t (fc)| 2 



This gives the Hamiltonian of Eq. [T] The quasiparticle 
energies E a (k) where a G {1,2,3} can be obtained by 
diagonalizing H^, and from the grand partition function 
Z = Tr exp[— /3?^mf], the grand potential ft \ 
is calculated as 



0(A) = 



-/9B a (fe) 



)• 



The ground state is self-consistently obtained by the 
order parameter A minimizing the grand potential f2, 
which in principle is equivalent to solving a gap equation 
but finds only stable solutions. 
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FIG. 1. The mixed-geometry lattice proposed and the phase diagram of exotic paired states at zero tempera- 
ture, a. Up-spin and down-spin components are selectively loaded on honeycomb (up-spin, sublattices A and B) and triangular 
(down-spin, sublattice A) lattices. The attractive Hubbard model is considered in the mixed-geometry lattice. Noninteracting 
band dispersions are sketched in b, where the gap in up-spin bands is tuned by the onsite energy modulation e (see Methods), 
c. The T-centered and K-centered Fermi surfaces are sketched. The high symmetry points T, K, and M are indicated in the 
first Brillouin zone. In a noninteracting system, a transition between the T- and K-centered shapes occurs at the van Hove 
singularity at = 1 where the Fermi surface is a hexagon (dotted line), d. The phase diagram is plotted as a function of 
up-spin density and polarization P — (n^ —n\)/ (n^ + n^) for the case of e = 0. The gapped and gap less paired phases and 
the normal phase are identified. The gap less states are further characterized by the Fermi surface topology denoted as z-FS-(X) 
indicating z Fermi surfaces being centered at the high symmetry points (X). Continuous (dotted lines) and discontinuous (solid 
lines) Lifshitz transitions are found, and the empty circle indicates the multicritical point. 
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FIG. 2. Characterization of paired states in the gapless phase. Representative Fermi surface configurations are 
shown with momentum-resolved order parameter over the first Brillouin zone (dashed line) in b, c, e, and f sampled from 
topologically distinct states given in a. The phase diagram (a) is plotted as a function of the chemical potentials ji^ and m for 
e = 0. The dotted and solid lines indicate continuous and discontinuous topological transitions which meet at the multicritical 
point indicated by the empty circle. The multicritical point corresponds to the fully filled down-spin band (/i^ = 0) in a 
noninteracting system showing a step-function-like singularity in the density of states. Momentum distributions in d and g 
further characterize the gapless paired states found in this system with 2-FS and 1-FS. The existence of plateaus of — 
in the shaded area enclosed by the Fermi surfaces is an indication of Luttinger's theorem. This mixed-geometry system shows 
pairing coexisting with polarization in momentum space in the gapless phase, in contrast to the phase separation in the usual 
breached-pair state. 




FIG. 3. Effects of the onsite energy modulation. The phase diagrams are examined for the onsite energy modulation 
e = — 1.0 (a), — 1.4 (b), and —1.8 (c). While the gapped phase is insensitive to e, the gapless phase (shaded) shows different 
Fermi surface topologies as e changes. Compared with the phase diagram for e = in Fig. [2^i, while the l-FS-(r) state is 
relatively robust, the areas of other incomplete-breached-pair states at e = are largely changed. In particular, a dominant 
phase at a large e is found to be the 2-FS-(r, K) state which turns out to be a different state from the one dominant at e = 0. 
This 2-FS-(r, K) state at finite e corresponds to case f in Fig. [4)3 rather than case d, showing a different momentum distribution 
with a fully polarized area as schematically shown in the inset of c. 
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FIG. 4. Schematics of pairing gap opening and Fermi surface formation. The band dispersions and spectral functions 
are shown along the K — T line in the first Brillouin zone. a. The noninteracting particle- like up-spin bands (<^) and the hole-like 
down-spin band (-£4,) are illustrated, possible points for level repulsion and gap formation are shown by empty boxes, b. The 
resulting quasiparticle bands with interaction turned on are sketched with possible locations of the zero excitation level (E — 0), 
showing the formation of zero (c), one (e), and two Fermi surfaces (d,f). The arrows indicate the possible evolutions from the 
noninteracting bands to the quasiparticle bands. The actual spectral functions of the down-spin component, A±(k, —oj), for all 
typical cases are selected for display in c-f. 



